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We consider bound states of heavy leptoquark-antiquark pairs (lepto-mesons) as well 
as leptoquark-antileptoquark pairs (leptoquarkonium). Unlike the situation for top quarks, 
leptoquarks (if they exist) may live long enough for these hadrons to form. We study 
the spectra and decay widths of these states in the context of a nonrelativistic potential 
model which matches the recently calculated two-loop QCD potential at short distances 
to a successful phenomenological quarkonium potential at intermediate distances. We also 
compute the expected number of events for these states at future colliders. 
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Particle physicists often lament the fact that the top quark is so short-lived. A long- 
lived top quark would have allowed us to study various aspects of the Standard Model 
(SM) that are otherwise inaccessible. For instance, had top quarks lived long enough to 
form bound states, they would have provided an exciting window into the nature of the 
strong interactions. But the SM top quark width, Tt->wb, grows as mf for large m t , 
and for mt = 175 GeV is about 1.55 GeV. This corresponds to a lifetime of approximately 
4x 10 -25 s — very short-lived indeed. However, there may still be a possibility of addressing 
these questions — not with top quarks, but with leptoquarks (LQs). These hypothetical 
particles carrying both lepton and baryon number appear naturally in many extensions 
of the SM. Some LQs may even be relatively light, with masses on the order of a few 
hundred GeV. (There has been a recent flurry of interest in these models due to data 
from HERA |]J which initially seemed to hint at the existence of a 200 GeV leptoquark, 
although new results from the Fermilab Tevatron || make this scenario unlikely.) In what 
follows, we will restrict ourselves to models in which each leptoquark has at least one 
trilinear coupling with some ordinary quark-lepton pair.0 Such leptoquarks will be color 
triplets having a spin of either zero or one. Their widths, which only grow linearly with 
their masses, will be proportional to the square of the above unknown coupling strengths. 
These couplings cannot be too large otherwise leptoquarks would have already been seen. 
Indeed, for reasonably large masses and small couplings, leptoquarks will have widths 
between about 1 and 100 MeV. As we will see below, this leaves plenty of room for the 
formation and observation of leptoquark bound states, either with another leptoquark 
(leptoquarkonium, LQ-LQ) or with a «, d, s, c or b quark (lepto-mesons, LQ-q)i It may 
even be reasonable to consider LQ-i states. 

The first step in identifying signals of these bound states is understanding their spec- 
tra. As with quarkonia, this is most easily done in a nonrelativistic potential model 
approach. Since leptoquarks are color triplets, the static QCD potential between an LQ 
and an LQ, or between an LQ and a q, is the same as that between a q and a q (although 
since leptoquarks have integral spin, relativistic corrections such as hyperfine splittings 
will be different). However, the region of the potential probed by our new bound states 

1 It is also usually assumed that a given leptoquark can only couple to quarks and leptons from 
a single generation — hence the terminology first generation leptoquark, etc. We instead make 
the less restrictive assumption that a given leptoquark couples to quarks from only one generation 
and leptons from only one generation, but the two need not be the same. 

2 Previous mention of bound states containing leptoquarks can be found in QQ. 
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will be somewhat different than that probed by the known quarkonium states. Thus, we 
cannot blindly use a generic phenomenological potential which is tailored only to the lat- 
ter, but must attempt to construct one more suited to our needs. Of course, this potential 
should be consistent with the known results from QCD as well as what we have learned 
from quarkonium studies. We will construct such a potential below (in the spirit of ||) 
and then use it to evaluate the energy splitting and the square of the wavefunction at the 
origin, |-0(O)| 2 , for the IS and 2S states of lepto quarkonium and lepto-mesons. This will 
allow us to calculate approximate widths and production cross sections, and estimate the 
expected number of events (on resonance) for these states. 

What do we know about the global structure of the static QCD potential V(r)1 At 
"short distances" it is quasi-Coulombic, with a perturbative expansion in the running 
coupling a s . The "long distance" behavior, describing quarks bound by tubes of chromo- 
electric flux, is expected to be linear. There is currently no known analytic method to 
describe V(r) between these regimes. However, many phenomenological potential models 
have been developed to fit the cc and bb spectra which basically probe this region, all of 
which have a quasi- logarithmic form at these distance scales. Finally, we note that V(r) 
is known to be concave downward for all r; that is, V'(r) > and V"(r) < §. 

We first discuss the short range potential Kihort ( r ) = — f-^rj where a v can be 
expanded in terms of a s . Until recently, this expansion was only known to one loop in the 
presence of rif massless flavors of quarks [jTj, although it was known to two loops in the 
absence of fermions [|8| . However, the full two loop potential has recently been computed 
||. The result, of course, is dependent on both the renormalization scale (fx = fi(r)) and 
scheme. We will use the BLM method [|10| for which 



Here the computation is performed in the MS scheme, and the scale at each order is chosen 
so that the coefficients Cj are n/-independent. Combining this with the general results of 
H yields c x = -8 and c 2 = 14/3 + 54tt 2 - 9tt 4 /4 - 220C 3 , as well as 

fi*(r) = -e I 1 6 6 V 4 - )\; fi**M = -e V 16 J (2) 
r r 

for the first two BLM scales (where 7 ~ 0.57722 and C3 ^ 1.20206). To two loops, the BLM 
method does not determine the scales fj,*** in (|IJ) and fx in (|2|). For simplicity, we choose 



both of these scales to be //**. An approximate value for a^g at any scale \x can be found 
from a reference value a^g(q) (we use cx-^s(m z ) = 0.118) and the following formula 
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which, like the potential, is complete to third order. Here 
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This form of the potential is only valid at energy scales far above quark pair creation 
thresholds, where all n/ quark flavors can be considered "light". Since this will be far 
from true in the applications that follow, we need a way to incorporate the effect of quark 
masses. A naive method would be to make rif a discontinuous function of /i, increasing by 
one unit as each threshold is crossed. However, this is not very realistic. In the absence of 
a full two loop calculation with massive quarks, we shall instead model these effects in the 
manner suggested in [jll[, letting n/ be an analytic function of fx: nf(fx) = V 



Here, the sum is over all quark flavors and the m g 's are the current quark masses (not to be 
confused with the constituent masses used later). As in [TJ] we choose m u = 0.004 GeV, 
m d = 0.008 GeV, m s = 0.2 GeV, m c = 1.5 GeV, m b = 4.5 GeV and m t = 175 GeV. (Our 
final results are rather insensitive to these numerical choices.) We believe this definition 
of V r s hort(^) to be valid out to r ~ 0.3 GeV -1 . 

The "intermediate distance" range is the least understood (from first principles) por- 
tion of the QCD potential. Running from the edge of the perturbative short distance 
region out to the long distance linear regime, it includes the quasi-logarithmic range 
(r ~ 0.9-5 GeV -1 ) where the low-lying cc and bb states live. Since the properties of cer- 
tain bound states containing leptoquarks will be sensitive to the transition between the 
perturbative and logarithmic regions, we want to connect them as smoothly as possible. 
To this end we would like an intermediate potential which attaches to our perturbative 
potential at some reasonable distance scale r a in such a way that V, V and V" are con- 
tinuous at r = r a . We also require that V(r) does a reasonable job of fitting the known 
(spin-averaged) cc (IS, 2S and IP) and bb (IS, 2S, 3S, IP and 2P) states. We choose to 
match at r a = 0.2 GeV -1 , which is safely (but not too far) within the perturbative region. 
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Since the first two derivatives of V will be continuous here, our results are somewhat in- 
sensitive to this choice. We choose the form of the intermediate potential to be a slightly 
modified version of the standard log potential [O: 



V iDt (r) = A\n 



To 
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As in |12[] we take A = 0.733 GeV. The constants ro, a and b are determined by our 



continuity criteria to be r ~ 5.501 GeV -1 , a ~ —0.013 GeV and b ~ 49.30 GeV. The 
second term in (]5|) is non-negligible only in the region r ~ 0.2-0.4 GeV -1 , so we do not 
expect it to have much of an effect on the quarkonium states. 
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Figure 1: Our version of the static QCD potential. 

As it stands, our potential does not have the appropriate linear behavior at "large" 
distances. We can, of course, fix this in any number of ways. However, we feel that this 
is unnecessary for our purposes since most of the states in which we are interested are not 
sensitive to this linear regime. The few states that we consider which have somewhat small 
reduced masses ( < 0.5 GeV) still lie mostly in the region well-described by the log potential 
(at least in the IS" case — see, for example, |13|]). In any event, the whole nonrelativistic 
approach of this paper is suspect for these states. We only include them in order to show 



certain qualitative trends.! Our full potential (Figure 1) is analytic except at r = r a , and 
is everywhere concave downward. 

Using the above potential and the nonrelativistic Schrodinger equation, we have nu- 
merically computed the bound state energies and wavefunctions at the origin for the IS 1 
and 2S states of lepto-mesons and leptoquarkonium (Table 1), for various constituent lep- 
toquark masses between 300 GeV and 1 TeV. (We also checked that V(r) fits the low-lying 
cc and bb states with approximately the same x 2 as the standard log potential alone.) 





System 


l^is(o)| 2 


I^2S(0)| 2 


AE2S-1S 


300 


LQ-u/d 


0.020 


0.010 


0.589 




LQ-s 


0.035 


0.018 


0.589 




LQ-c 


0.190 


0.097 


0.589 




LQ-b 


1.285 


0.649 


0.593 




LQ-i 


1,985 


308.5 


1.346 




LQ-LQ 


4,574 


640.8 


1.751 


500 


LQ-c 


0.191 


0.097 


0.589 




LQ-b 


1.299 


0.656 


0.593 




LQ-i 


3,072 


445.8 


1.542 




LQ-LQ 


1.819xl0 4 


2,622 


2.725 


700 


LQ-c 


0.191 


0.097 


0.589 




LQ-b 


1.305 


0.659 


0.593 




LQ-i 


3,791 


538.1 


1.649 




LQ-LQ 


4.489 xlO 4 


6,776 


3.604 


1000 


LQ-c 


0.191 


0.097 


0.589 




LQ-b 


1.310 


0.662 


0.593 




LQ-i 


4,486 


629.2 


1.741 




LQ-LQ 


1.167xl0 5 


1.820xl0 4 


4.808 



Table 1: Values of \ip n s(Q)\ 2 and the 1S-2S energy splitting for 
lepto-mesons and leptoquarkonium. All quantities are in appro- 
priate GeV units. 

This analysis also required choices for the constituent quark masses. For the c and b 

1.5 GeV and m b = 4.906 GeV. For the top 



quarks we used those from |12 



3 We have also ignored the effects of leptoquark-Higgs interactions in our potential. These very 
model-dependent couplings, even if large, would only affect the LQ-LQ (or possibly LQ-i) states 
- and even then only if the Higgs mass is "small enough". We neglect the model-dependent 
leptoquark couplings to electroweak gauge bosons as well. 
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quark we used rat = 175 GeV, while for the lighter quarks we chose m s = 0.5 GeV 
and m u = rrid = 0.35 GeV. The states LQ-u/d and LQ-s are only presented for a lepto- 
quark mass of 300 GeV, since the reduced mass of the system is basically independent of 
the leptoquark mass in the range of interest. The radii of the IS" heavy-heavy states in 
Table 1 (LQ-LQ and LQ-i) lie well within the perturbative regime of V(r), ranging from 
approximately 0.02 to 0.05 GeV -1 as the reduced mass decreases. The 25" heavy-heavy 
states have radii from about 0.05 to 0.22 GeV -1 . Here we begin to see some sensitivity to 
the short-intermediate matching region for lighter leptoquark masses. The IS LQ-b states 
all lie around 0.6 GeV -1 , and are also somewhat sensitive to the matching. All other states 
considered have radii of 1.2 GeV -1 and above, and thus mainly see the logarithmic part 
of the potential. (For an alternative treatment of LQ-c and LQ-b states, see [|].) 

Technically, the results in Table 1 represent spin- averaged quantities since we have not 
included hyperfine effects in our potential. However, for scalar leptoquarks this statement 
is trivial since their leptoquarkonium states can only have spin-0 and their lepto-meson 
states spin-1/2. For vector leptoquarks there are hyperfine splittings. These are negligible 
for the lepto-mesons containing lighter quarks, but larger for the LQ-i states and for 
leptoquarkonium. However, they are still substantially smaller than the 1S-2S splittings 
( < 0.1 Ai?2S-is), justifying our assumption that the static potential alone is a good 
starting point. For simplicity, in the remainder of this paper we will only consider scalar 
leptoquarks. 

Our analysis so far has only utilized the mass and spin of the leptoquark, along with 
the assumption that it transforms as an SU(3) C triplet. The study of leptoquark bound 
state production and decay, however, is more model-dependent; for example, 77 — > LQ- 
LQ involves the leptoquark's electric charge. We follow [|14| in assuming that at the mass 
scales considered, the SM is extended only by the addition of leptoquarks which respect 
baryon and lepton number and SU(3) C x SU(2)l X U(1)y- The leptoquark isospin and 
hypercharge then determine its coupling to the SM gauge bosons, and to which quark- 
lepton pairs it may couple. We explicitly consider only one (presumably the lightest) 
scalar leptoquark. A degenerate weak isospin doublet or triplet of leptoquarks would not 
drastically change our conclusions, but could enrich the phenomenology (by giving rise to 
the decay LQ-LQ — > W + W~ , for example). Except in the case of £ + £~ — > LQ-LQ, we 
assume a purely chiral leptoquark coupling; if its trilinear coupling to a lepton of chirality 
i (i = L,R) and quark of opposite chirality is denoted by A^, then we assume that either 
Xl 7^ 0, Xr = or Xr 7^ 0, Xl = in order to satisfy magnetic moment constraints (see 
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below). For convenience, we write Aj = eAj, with e 2 /An = a em {m z ) = 1/128. To ensure 
that the bound state widths are substantially smaller than the 1S-2S splitting (making 
the states easily identifiable), we require weakly coupled leptoquarks (Af < 1). This leads 
to bound state widths of at most (9(100 MeV). Finally, we continue to ignore the model- 
dependent leptoquark couplings to W, Z and Higgs bosons - - including the effects of 



partial widths such as LQ-LQ — > ZZ would not greatly change our results. 

The production cross-section (on resonance) of a narrow width bound state B from 
an initial state i, summed over all decay modes of the resonance ||15|| . is well-approximated 
by o"t t = (lQiiS/Mg) ( r ^T' ). Here, V and Ts^i are the total and partial widths of the 
bound state, and the rest-frame total energy E is well within V of the resonance mass 
Mb- In the particular case of unpolarized incoming beams, the spin- multiplicity factor 
S is given by Nb/ (N1N2), with Nb the spin-multiplicity of the resonance and Ni,N% 
those of the incoming particles comprising the state i. When the colliding beams are not 
monochromatic compared to T, we assume a gaussian distribution in E strongly peaked 
at the resonance mass Mb with a spread of a E (a E > T). At the energies considered, 
a tot = (l67ri?/Mg) (a/1F) y ^E j' ^ e assume that leptoquarks would first be discovered 
in non-bound state production, and that m LQ would be sufficiently well-determined to 
allow a search for the resonances considered here. 

We first treat lepto-mesons {Mb — rn LQ + m q ). Throughout, our notation will not 
distinguish between particle and antiparticle so that (for example) LQ-q represents LQ-q 
or LQ-q, and TLQ^ eq represents Y L q^ e ± q or Tj^^ e± -. The principal lepto-meson decay 
channels should be the spectator decays of the quark (J? B^>(LQ)q +x = ^q^q'+x) an d 
leptoquark: 

Y B ^iqq = T LQ ^e q = aem ^ 1 I — — li J mLQ (£ = e ,^,r). (6) 

4 \ m LQ J 

(We will ignore leptoquark couplings to neutrinos as well as all lepton masses.) As men- 
tioned earlier, for large m LQ the leptoquark decay width increases only as m LQ . The 
quark decay width T B^(LQ)q +x * s negligible except in the case of the top quark, where 
r B ^{LQ)q +x = Ft->wb- We will also need the very narrow width FB^£ alb = ^ + 4 °^ rg^ 7 
for decay into £7 with helicities a/2 and b respectively: 

r B ^ = * (4 TOe ,„A,) 2 \mf {%- - — ■ m 




16n V / \M B m q J m LQ 



Here N c = 3 and Q is the electric charge (in units of the proton charge). The width 
Ts^i-y is inconsequential as a contribution to the total width, which is well-approximated 
by T = r * B^(LQ)q '+x + ^B^iqq- Over the mass range we consider, the LQ-t width Tb^zu 
ranges from A| ■ 255 MeV to Af • 1835 MeV, while for lighter quarks the LQ-q widths 
Y B^iqq all range from about Af ■ 586 MeV to Af ■ 1953 MeV. We may now consider the 
production of lepto-mesons at an e~7 collider (for non-bound state production, see Hl6"l ). 
With an expected beam spread of a E = 0.05-s/s ~ 0.05-Mb [0, the cross-section for 
a polarized e~7 collider is cr^ i ab = (16tv/M b ) (^/f) The spin-multiplicity 

factor S has been specified implicitly in this case (S = 1) as is done below for the other 
cross-sections computed. For unpolarized incoming beams, S = 1/2. 

With the exception of LQ-t, for which the top decay width dominates, the 1S-2S 
splitting is larger than the width of the states for Af < 1/10. (For LQ-t, there may still 
be the possibility of resolving the 1S-2S peaks for very heavy leptoquarks — see Table 1.) 
In Table 2 we present the expected number of events per Af for the production of IS 
states. Given a yearly integrated luminosity of up to 50 fb _1 



17 , these resonances could 



be observable. Resolving these states, however, will be difficult in practice given the broad 
collider energy distribution and the non-trivial task of precisely reconstructing the lepto- 
meson given that there is a jet in the final state. 





LQ-t 


LQ-b 


LQ-c 


LQ-s 


LQ-u 


LQ-d 


300 


72.36 


13.58 


98.51 


40.36 


194.5 


48.11 


500 


26.63 


3.127 


21.31 


8.765 


41.98 


10.43 


700 


12.24 


1.171 


7.769 


3.202 


15.30 


3.807 


1000 


4.979 


0.410 


2.665 


1.100 


5.246 


1.307 



Table 2: Estimates of the number of events/A? for IS" lepto- 
mesons produced approximately on resonance at a polarized e~7 
collider. Masses are in GeV. We assume a E = 0.05-Mb and an 
integrated luminosity of 10 fb _1 . The particle-antiparticle assign- 
ments shown in the column headings were chosen to maximize 
production. 



the partial width for the spectator decay 
= 2r^Q^^ 9 . For a very weak trilinear 



Turning to leptoquarkonium {Mb — 2m LQ ), 
of either constituent leptoquark is F B ^eq(LQ) 
coupling (A^ < 1/1000) this decay channel is not dominant. The two gluon or two photon 
(with helicities a and b) decay modes are more significant: 



B^g a 9b 



Oab 

2nN r 



Ml 



(8) 



_ _ B 



On the other hand, for A? > 1/100 we have V ~ r s _^ g (£Q). We can now discuss the 
production of leptoquarkonium at a 77 collider. Assuming a 77 center of mass energy 



resolution a E of 0.15^ ~ 0.15Mb [[1711 , the cross-section for unpolarized incoming photons 



is 



atot= M|V8^r- (10) 



(There is an additional factor of 2 in this equation to compensate for the factor of 1/2 
that was included in Eq. (|]) to account for the identical photons in the final state.) 
Eq. (^) shows that production will be greatly enhanced or suppressed depending on the 
value of Qlq (\Qlq\ = 1/3, 2/3, 4/3 or 5/3). We have taken the most optimistic possibility, 
IQlqI = 5/3, in Table 3. It is interesting to consider a very small A?, say Af = 1/10000, 
where the dominant decay width (~ 1.7 MeV for m LQ = 300 GeV) is into gluon pairs. In 
this case the 77 decay channel, with a partial width of 0.45 MeV, has a branching ratio of 
about 1/5. The resulting narrow resonance in 77 — > 77 will stand out in a spectacular way 
over the continuum background [IB], which should be more than an order of magnitude 
smaller (at somewhat less than 1 fb within a reasonable mass bin). Similar results are 
obtained even for |Qz,q| as small as 2/3. 

Finally, for leptoquarkonium production at an e + e~ or (J. + n~ collider we will need the 
decay width into a pair of charged leptons (i + i~) with helicities a/2 and 6/2: 

/ ~ ~ \2 |^(0)| 2 m2 
r B ^+ r = W27rAf c (a em A +a A_ a ) /J1/f2 i a^.w (H) 



(For a review of non-bound state production, see [|n§.) The cross-section for a polarized 
iH- coUider operating at the _ee is then , tot , o6 = (16./MD (^i), which 
vanishes unless we relax our condition that the coupling be purely chiral; in fact, it is 
maximized if we take A^ = Xr (as in Table 3). The maximal cross-section will be for an 
LQ that couples to it, which we assume here. Since the next generation e + e~ colliders are 
expected to attain at best a E ~ COly^ 3> T, and the measured electron magnetic moment 
requires an extremely small value for Xl^r [|0[ (and thus a tiny partial width r B ^ e + e -), 
production at these machines would be heavily suppressed. By contrast, proposed muon 
colliders could attain a E ~ 2.1 x 10 _5 -y/i [^TJ. Furthermore, the constraints from the muon 
g — 2 are much looser, allowing Af, = \ 2 R as large as 0.14 for m LQ > 300 GeV. Results 
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for — > LQ-LQ may be found in Table 3. In this case, the useful decay mode would 
most likely be B — > \i~ q^ + q (with a branching ratio of 0.97 for m LQ = 300 GeV) so that 
this channel would not only be visible, but might also have the potential of resolving the 
15-25 splitting. 





H + H — > LQ-LQ 


77 -> LQ-LQ 


300 


578.0 


853.8 


500 


52.55 


264.0 


700 


12.05 


121.1 


1000 


2.568 


52.91 



Table 3: Estimates of the number of events (on resonance) at 
polarized fx + \T~ and unpolarized 77 colliders for 15 leptoquarko- 
nium. Masses are in GeV and the integrated luminosity is 10 fb _1 . 
In the fi + \i~ case, we assume that a E < V, \Qlq\ = 1/3, and LQ 
couples to \it with Af, = X% = 0.1. In the 77 case, we take 
a E = 0.15M B and \Qlq\ =5/3. 

In conclusion, we have constructed a version of the static QCD potential which 
matches (a slightly updated form of) the recently calculated two-loop perturbative piece 
at short distances to a successful phenomenological quarkonium potential at intermediate 
distance scales. We then used this potential to determine the spectra and decay widths of 
bound states containing heavy leptoquarks; namely, lepto-mesons and lepto quarkonium. 
The expected number of events for these states at future colliders were also estimated. 
For reasonable values of leptoquark masses and couplings, many of these exotic hadrons 
should be observable. 



Note: After the completion of this work, we noticed that a similar version of the static 



QCD potential has recently been constructed in [[22] , although used for different purposes. 
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